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INTRODUCTION 


The purpose of this study was to develop a model for worldwide 
cloud cover using a satellite data set containing infrared radiation 
measurements. Other cloud models exist [1-3]. These early cloud 
models used primarily ground-based cloud observations. The satellite 
data set containing Day IR, Night IR, Incoming, and Absorbed solar 
radiation measurements on a 2.5-degree latitutde-longitude grid 
covering a 45 month period of record has recently become available. 
There was originally a 2-year period of similar data on an NMC grid. 
The first step was to convert these infrared data to estimates of 
cloud cover. The statistical analysis of classification of cloud 
region characteristics was then performed. 

There are several reasons for desiring a cloud model based on 
satellite data. The ground-based data are much more limited in scope. 
Some fairly large areas of the world have either no data or very 
sparse data, and models using ground-based observations necessitate 
a number of assumptions, including on occasion that a region is 
essentially like its antipodal location. A good worldwide cloud 
cover model is needed for the purpose of studying the relationship 
between cloudiness, precipitation, and the Earth radiation budget. 
CONVERSION OF SATELLITE IR MEASUREMENTS TO CLOUD COVER 

A major initial task was to derive cloud cover estimates from 
the satellite infrared data. The method used in this investigation 
follows the suggestions obtained through personal communications with 
Thomas I. Gray, Jr. |4j. 

Albedo is defined as the reflective power, or the fraction of 
incident light that is reflected by a surface or body. Included in 
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the satellite data are the amount of incoming solar radiation, and 

the amount of absorbed solar radiation, I . . The satellite observed 

an 


albedo A, is estimated by 


A = - I 

in an in 


( 1 ) 


If the Earth's surface absorbed all solar radiation, then the cloud 
cover might be taken simply as 1 minus albedo (assuming also that 
clouds reflect all solar radiation). Different parts of the Earth's 
surface, however, have differing radiances. For example, the albedo 
of the ocean is approximately 5 percent (95 percent of the solar 
radiation being absorbed), while tin* Sahara desert reflects approxi- 
mately AO percent of the solar radial ii»n reaching it. 

To determine cloud cover, we iweded lo obtain the background 
radiation of the region of the Earth of inlerest. To do this, for a 
given season and a specified location, we calculated A from equation 

(1) for every day of a season and observed the minimum value, A . • 

min 

This minimum value should occur on I lie day of least (hopefully near 

zero) cloud cover. If r is the rel l»'cl.ince of Ihe clouds and x is 

the fraction of cloud cover, then the basic formula may he written as 

A = X r » (l-x) * A . 

min 

from which we have the fraction ol cloud (over, x, as: 


X = 


A - A . 

mm . 

' ^miii 


( 2 ) 


This forsHila requires a knowledge ol r which varies. 

A way to estimate the cloud rel Irctaiice r, is by observing the 
difference between the Karth's surface temperature and the tempers* 
ture equivalent of the satoll ite-ohserved day time infrared reading 
(denoted by IRq)« The radiance of the IR|^ by Stefan's law is equal 
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to 5.75 X 10 T (watts/m ), where T is the temperature equivalent 
in degrees Kelvin. Putting z = (surface temperature - T), (units 
degree Kelvin), the following relationship has been observed: 

r = -0.000265z^ + 0.0295z + 0.10. (3) 

A surface temperature of 30° C for latitudes within 25 degrees of the 
equator and -5° C for latitudes within 25 degrees of the pole was 
used. Interpolations were used for intermediate latitudes. 

Using this method, the proportion of cloud cover was calculated 
for each day of the year over the 4-year period covering the entire 
globe. We considered the data for the four seasons separately, and 
we developed a separate cloud cover model for Winter, Spring, Sumner, 
and Fall, where Winter consists of the months December, January, 
and February, etc. 

Because this investigation is based on derived cloud cover 
estimates and may be subject to criticism, it is noted that ground- 
based cloud observations are also estimates as well as cloud cover 
obtained by satellite photography. We make this conjecture: Those 

variables which are not well defined in the IR to cloud cover 
conversion procedure will have small contributions to climatic model- 
ling of the clouds over the entire season. For a specific day and 
area the preceding procedure may not be entirely satisfactory for 
synoptic cloud cover analysis. It should be noted that because of 
the orbit of the TIROS satellite the estimates in the near polar 
regions are degraded somewhat. 

A PROBABALISTIC MODEL FOR CIX)U0 COVER 

The proportion of cloud cover over any grid square is a random 
variable which has some probability distribution associated with it. 
Falls (3) found that the Beta distribution could be used to represent 
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the probability distribution of the proportion of cloud cover. 
Hend?rson-Sellers [5] has also found the Beta distribution useful as 
a model for the probability distribution of the proportion of cloud 
cover. The Beta probability density function with parameters a and 
b is given by: 


f(x) 


r(a + J>)_ * (1 - x)** ^ 

F(a)T(b) 


for O^xSl, a>0, b>0 

The Beta probability density function can assume a variety of 
shapes. It can be mound shaped, U-shaped, or J-shaped with varying 
amounts of skewness. Table 1 shows tlie relationship between the a 
and b parameters and the shape of the freipiency curve. 



TABLE 1: Shapes of the Beta Urohability Density Function 

For Different a and b Parameters 

Shape Parameters 

Hound a > 1 b > 1 

J a > 1 b < I 

Reverse J a < 1 b > 1 

U a < 1 b < 1 

Uniform (0,1) a = 1 b = 1 

The question arises as to how well the Beta distribution actually 
fits the cloud cover data at hand. To answer this question we first 
estimated the a and b parameters using the method of moments. These 
estimates are given by 


b = (1 - x) |x (1 - x) - s^)/s^ 
a = X b /{I - x) 

m 2 

where x is the sample mean proportion of cloud cover, and s is the 
sample variance of the proportion of cloud cover. Next, we 
constructed histograms for the cloud cover at several grid locations. 
A grid location was selected at random from latitude circles 12.5^ 


! 

|. 


i - 
i ■ 


f 





apart beginning with 50^ N latitude and extending to 62.5^ S latitude. 
A histogram was constructed for each selected location on the basis 
of the cloud cover data for the Winter quarter (December, January, 
February) for the four years of data. Each of the histograms were 
constructed on the basis of approximately 350 cloud cover values. 

The corresponding Beta curves are shown superimpt>sed on the histo- 
grams in Figures 1-2. 
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A WORLD WIDE CLOUD COVER MODEL LSIM', THE RETA DISTRIBUTION 

The (a,b) paramclers of the Beta di;:tribut ion f?ive a good deal 
of information about the cloud cover characteristics of a given 
location as illustrated in the previous table and figures. Thus, 
we used the estimated parameters (i,b) in determining regions of 
homogeneous cloud cover. 

As a first step, it ia instructive to consider the frequency 
histograms of the calculated (a,b) values for each of the four 
seasons. There arc i0,224 grid points over the globe, and each grid 
point has an (a,h) pair associated with it lor each season. The 
following frequency histogr.ims (Kiguies l-<») of tlu* 10,224 (a,b) 
pairs give some indicat i<>n ol where the par.imeter pairs are falling 
in the (a,b) plane. 
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Figure 3: Vinter Frequency llislogr:im for Global (a,b) Values 
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Figure 6: F.ill Frrqunuy Hislfgram lor Uloh.il (.i,b) Values 
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The numbers in the above I roqin-m y bislograms indie.ite how many 

A. A 

(a,b) paraiaelers over the entire glohr tall in the specif ie’ Mock. 

A A 

For example, there are 644 (a,b) pairs ^lutially siieti that 

A A 

1.0 < a ^ 1.5 and b > 5.0 in the Winirr qiMiler. 

The 36 regions in the (a,b) pl.jiio (row the fre<|uency histogram 

A A 

for global (a,b) values form a basts lor <i^aermining homogeneous 

A A 

cloud cover regions. Grid points «>ii the globe which have (a,b) 
parameters falling in tlie saau; blot k have very similar cloud cover 
characteristics. This concept leads to a preliminary cloud cover 
swdel. Each grid point on the globe is assigned to one of 36 groups 

A ^ 

depending on the bl'jck in which (a,h) falls. A FORTRAN program was 
used to label each grid point with 0 • 9 or A ' Z depending on the 
group the grid point fell into. These labels were printed in a 
rectangular array maintaining the latitude and longitude position of 
each point. This procedure results in a m.ip o. the globe containing 
a large number of contiguous regions \diich have thf :iar..e nasic cloud 
cover characteristics. The resulting maps are si^nwii in Figures 7*10 
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Figure 8: Global Cloud Cover Classification 

Based on (a,b) Parameters of the Beta 
Distribution - Fall 
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with global maps superimposed. These maps have been simplified by 
combining several of Llie 36 blocks and recoding the maps. The 
recoded maps are shown in Figures 15-18. The key for the original 
and maps is given in Figure 11. The key for the recoded maps is 
given in Figure 12. 
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Figure 

Tne maps may be interpreted acconling to the distributional charac- 
teristics of the various regions. Typical Beta frequency curves for 
the 9 recoded regions are given in Figures 13-lA. 
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Figure 12: Key For Recoded Homogeneous Cloud Cover Regions Map 
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Figure 15: Global Cloud Cover Classification Based o 

(a,b) Parameters of the Beta Distribution 
“ Summer 






Figure 16; Global Cloud Cover Classification Based on 
(a,b) Parameters of the Beta Distribution 
- Fall 












Figure 18; Global Cloud Cover Classification Based on 
(a,b) Parameters of the Beta Distribution 
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WORLD WIDE CLOUD COVER MODEL USING MEAN AND STANDARD DEVIATION CONTOUR 
LINES 

One may wish to determine mean and standard deviation of the 
proportion of cloud cover at any given location. To do this we have 
developed contour maps which give the mean or standard deviation. 

The methodology is similar to that used to determine homogeneous 
cloud cover regions. We first calculated the mean and standard 
deviation of cloud cover for each grid point for a particular season. 
To obtain a contour map for the mean cloud cover for a particular 
season we used s FORTRAN program to print the first digit of the 
mean cloud cover for each grid point on the globe maintaining the 
latitude and longitude position. Contour lines were then traced on 
a transparent map overlaying the printout. Separate maps were 
constructed in the same way for the standard deviation of cloud cover. 
These maps are shown in Figures 19-26. 

The maps in Figures 19-26 may be used as a general guide to 
the cloud cover characteristics for any place on the globe. Also, 
the aiean and standard deviation of the proportion of cloud cover 
obtained from these maps may be used (in the absence of the pre’*iously 
given auips in Figures 7-10) to obtain estisuites for the parameters 
(a,b) in the Beta distribution. 

CONCLUSIONS AND RECOMMENDATIONS 

The current cloud cover model illustrates a useful objective 
methodology for cloud cover classification. The developed maps can 
be useful for those who need sosm inforaution regarding the cloud 
cover characteristics for any particular location on the globe. 

This is particularly useful in that the practicing clisMtologist can 
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Figure 20: Standard Deviation Cloud Cover - Summer 

Printed x 10 => Standard Deviation 







Figure 21: Mean Cloud Cover - Fall 















Figure 23: Mean Cloud Cover - Winter 




} 


IK 


t 



u 

0 ) 

> 

o 

a 


o 


u 


VM 

o 

0 

o 

•H 

4 J 

CT3 

•H 

> 

<D 

O 

T 5 

U 

T 3 

C 

cd 


t/3 


-sf 

CNJ 


<D 

t-l 

D 

60 


Pm 


29 


Winter 








t 

. ^ 


Spring 




32 


obtain a great deal of cloud rover information without going through 
large volumes of data. 


There are problems with the curreMit cloud cover model. The 
most obvious of these is the lack of data. To make a good climato- 
logical model a reasonably long record length is required. The 
satellite data available for this study comprise approximately 44 
months. This is sufficient lor some model development, but a longer 
period of record would be desirable. Also, temporal persistence 
cuts down on the actual number of independent cloud cover observa- 
tions. Another problem is that the cloud cover used in the model 
is derived. This is not to say thal the cloiul cover values used 
are not accurate, but the derived cloud cover should be compared 
with some independently observed cloud cover. It should be noted 
that any type of cloud cover measurement will have some degree of 
error associated with it. 

The topic of error analysis in the estimation of the parameters 
(a,b) was not covered for several reasons. There are several types 
of errors which made the problem quite complex. First, as previously 
mentioned, there is an unknown error component in the derived cloud 
cover itself which may vary from place to place. Also, there is 
error in recorded satellite immsuremenls. The cloud cover values 
have a temporal correlation which is not fully known. The short 
length of the record certainly contriliules to the error (the degree 
of error depends on the temporal persistence of cloud cover and the 
long term cyclic behavior of cloud cover). Another factor is that 
the satellite measurements were taken for a given grid location at 
the same time each day. 



It is worth remarking that even with the above*mentioned errors, 
it is still possible to use the data at hand to classify areas of the 
globe into homogeneous cloud cover regions. The methodology is sound, 
and the presently developed models should give a good overall picture 
of global cloud cover characteristics. 

The above comments lead to three major recommendations. First, 
the derived cloud cover should be verified using sosw independent 
measurements. Second, spatial and temporal persistence of cloud 
cover need to be investigated. Third, the model should be updated 
with a longer length record. 
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